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Abstract 


Severe Acute and Middle East Respiratory syndrome coronaviruses (SARS-CoV and MERS- 
CoV) encode multifunctional papain-like proteases (PLPs) that have the ability to process the 
viral polyprotein to facilitate RNA replication as well as antagonize the host innate-immune 
response. The latter function involves reversing post-translational modification of cellular 
proteins conjugated with either ubiquitin (Ub) or Ub-like interferon stimulated gene product 15 
(ISG15). Ubiquitin is known to be highly conserved among eukaryotes but surprisingly ISG15 is 
highly divergent among animals. The ramifications of this sequence divergence to recognition of 
ISG15 by coronaviral papain-like protease at the structural and biochemical levels are poorly 
understood. Therefore, the activity of PLPs from SARS-CoV, MERS-CoV and mouse hepatitis 
virus (MHV) was evaluated against seven ISG15s originating from an assortment of animal 
species susceptible, and not, to certain coronavirus infections. Excitingly, our kinetic, 
thennodynamic and structural analysis revealed an array of different preferences among PLPs. 
Included in these studies is the first insight into a coronoavirus PLP’s interface with ISG15 via 
SARS-CoV PLP in complex with the principle binding domain of human and mouse ISG15s. 
The first X-ray structure of the full-length mouse ISG15 protein is also reported and highlights a 
unique, twisted-hinge region of ISG15 that is not conserved in human ISG15 suggesting a 
potential role in differential recognition. Taken together, this new infonnation provides a 
structural and biochemical understanding of the distinct specificities amongst coronavirus PLPs 
observed and addresses a critical gap of how PLPs can interact with ISG15s from a wide variety 
of species. 

Ub, Ubiquitin; ISG15, Interferon stimulated gene product 15; PLpro, Papain-like proteases; PLP, Papain-like proteases (for viruses that have 
more than one PLpro); MHV, Mouse hepatitis virus; CoV, Coronavirus; vOTU, Viral ovarian tumor domain proteases; IFN, Type-I interferon; 
Ubl, Ubiquitin-like domain; NS and Nsps, Nonstructural proteins; USP, Ubiquitin-specific proteases; ITC, isothermal titration calorimetry 
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INTRODUCTION 


Coronaviruses (CoVs) are enveloped, positive-stranded RNA viruses that cause mild to 
severe infections in a wide range of mammals and birds. Specifically, severe acute and Middle 
East respiratory syndrome coronaviruses (SARS-CoV; MERS-CoV) are well-recognized viral 
pathogens that have emerged from different animal reservoirs to cause deadly disease in humans. 
SARS-CoV first emerged in 2002 with a case fatality rate of 10%, claiming the lives of over 800 
people and infecting more than 8000 ’ . Ten years later, MERS-CoV emerged with a shocking 
case fatality rate near 35% and has spread to 27 different countries to date . The continuing 
threat of MERS-CoV was recently underscored by one of its most recent outbreaks in Republic 
of Korea. This outbreak quickly led to 36 deaths, which has brought the total MERS-CoV global 
deaths to over 600 4 ' 5 . 

Similar to other positive stranded RNA viruses, CoVs encode two types of cysteine 
proteases, including the papain-like protease (PLP) and 3C-like protease, also known as the main 
protease. Together, these enzymes cleave the viral polyprotein into 16 different nonstructural 
proteins (Nsps 1-16) in order to generate the membrane-bound replicase complex for RNA 
replication 6 ' 7 ' 8 ' 9 . CoVs encode either a single PLP, termed PLpro, or two PLPs that process a 
total of three cleavage sites within the polyprotein ’ . For instance, SARS-CoV and MERS- 
CoV encode a single PLpro, while other CoVs such as mouse hepatitis virus (MHV) encode for 
both the papain-like protease 1 (PLP1) and papain-like protease 2 (PLP2) (Figure la). 

Beyond cleaving the viral polyprotein, PLPs have additional activities that promote virus 
replication. The X-ray structure of the first CoV PLP determined from SARS revealed that these 
enzymes resemble the structure of human ubiquitin-specific proteases (USPs) and are thereby 


4 




known as viral USPs, often acting as deubiquitinating enzymes with the ability to remove the 

o 

post-translational modification ubiquitin (Ub) from target proteins . Some PLPs are also 
delSGylating enzymes with the ability to reverse the post-translational modification of the Ub- 
like protein interferon stimulated gene product 15 (ISG15) from cellular proteins 1 . Such 
activities were implicated in SARS-CoV’s suppression of the innate immune responses, 
particularly antagonizing type-I interferon (IFN) signaling and chemokine and cytokine 
production ’ . Simultaneous disruption of both activities, either by mutation in MERS CoV 

PLpro’s Ub/Ub-like binding region 14 , or a destabilizing mutation in MHV PLP’s ubiquitin-like 
domain (Ubl) 1(1 has been observed to prevent antagonization of IFN signaling, chemokine and 
cytokine production, as well as viral pathogenesis respectively. Overall, these studies suggest 
that the deubiquitinating and delSGylating activities of PLPs, sometimes packaged together as 
deubiquitinating enzyme (DUB) activity, are a likely contributor to viral pathogenesis. 

Intriguingly, the precise role of each individual activity in promoting PLPs ability to act 
as an IFN antagonist has yet to be precisely defined. So far, several X-ray structures of PLpro 
and PLP2 bound with Ub molecules have been reported 8 ’ 15 ' 16 . However, no structure of any 
CoV PLP in complex with an ISG15 molecule has yet been reported. Due to the lack of 
structural detail about the interface between CoV PLPs and ISG15, understanding how CoVs 
PLPs specifically engage with ISG15s versus Ub has been especially difficult to understand. The 
lack of structural information has led to problems in defining the differences between CoV PLPs’ 
deubiquitinating and delSGylating activities among different CoVs. A further complication in 
our understanding stems from the fact that the sequence identity of ISG15 among mammals is 
low, ranging from 58% to as low as 35% when a broader range of animals is compared (Figure 
lb). In contrast, Ub is highly conserved among eukaryotes. 
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The sequence diversity of ISG15 among animals and the impact of this diversity on the 


recognition of these ISG15s by viral USPs has not been well studied despite the potential 
implications. For example, certain CoVs, such as SARS-CoV and MERS-CoV, are known to 
replicate and survive in a broad range of animals including bats, camels, mice, civets, shrews, 
badgers, pigs and humans ’ ’ . In contrast, other CoVs such as MHV only replicate in mice, 
suggesting that some CoV PLPs may have evolved strict specificity for their single host’s 
ISG15 19 . The potential for species-species variances in ISG15 have already been shown to be a 
key factor in other viruses including influenza B where it was shown that NS1 is unable to bind 
mouse ISG15 unlike its non-human primate and human counterparts limiting influenza B virus 
infection in mice ’ ’ . Also, biodiversity that occurs between species within ISG15s was 

recently shown to impact the ability of nairovirus viral ovarian tumor domain proteases (vOTUs) 
to effectively process certain ISG15s “ . Overall, many CoV PLPs have been observed to show 
robust deubiquitinating activities; n ’ 24 however, recent in vivo studies have started to unveil the 
importance of PLPs delSGylating activity 25,26 . This has led to a need for a better and more 
detailed understanding of the interactions between CoV PLPs and ISG15 at the molecular level. 

Adding to the ambiguity of how CoV PLPs or other delSGylating enzymes interact with 
ISG15, especially when it comes to understanding species-to-species variations, is the lack of 
available full-length ISG15 structures that have been determined. This is especially important 
when the intramolecular arrangement of the two domains of ISG15 is considered. For example, 
the hinge region of ISG15 shows a significant degree of sequence diversity between the different 
species (Figure lb). Although several structures of human ISG15 (hISG15) have been reported, 
no complete structure of ISG15 from another species has been resolved. Moreover, the recent 
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elucidation of the structure of the C-terminal domain portion of mouse ISG15 (mISG15) 

'j'l 

prompted questions on the impact of ISG15 biodiversity on the overall ISG15 structural fold . 


To address these critical gaps in our understanding of ISG15 recognition by CoV PLPs, 
PLpro from SARS-CoV and MERS-CoV, as well as PLP2 from MHV, were examined for 
differences in their selectivity among ISG15s from various animals. The X-ray crystal structures 
of SARS-CoV PLpro bound to the C-terminal domain of ISG15 originating from mouse and 
human were determined and analyzed in conjunction with enzyme kinetic and thermodynamic 
data derived from ITC. In addition, the first X-ray structure of a complete non-human ISG15 
structure, mouse ISG15 (mISG15), was also determined to elucidate potential sequence and 
structural differences that may account for species specificity of CoV PLPs. Together, these 
studies provide significant and new insight into the CoV PLP’s ability to accommodate the 
structural differences not only between Ub and hISG15, but also different species ISG15s. 

RESULTS 

Species specific cleavage of proISG15s by CoV PLPs. 

Currently, there is a paucity of available biochemical data on the impact of species-to- 
species sequence variations within ISG15 and the effect that these differences may have on the 
ability of CoV PLPs to recognize and cleave ISG15. Therefore, we employed a recently 
developed assay for delSGylating vOTUs that takes advantage of the ability of a protease or 
DUB to cleave immature ISG15 23 . The PLpro enzymes from SARS-CoV and MERS-CoV, as 
well as the enzyme PLP2 from MHV, were purified and assessed for their ability to cleave 
proISG15s derived from 7 different species including human, sheep, northern tree shrew, 
jackknife fish, mouse, dromedary camel, and vesper bat (Ligure 2). Each of these proISG15 
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proteins is appended with the proISG15 extension from Homo sapiens (Figure 2). ISG15 from 


jackknife fish was included in the analysis to add a more distantly related ISG15 homologue. 

SARS-CoV and MERS-CoV PLpros are both capable of fully processing proISG15 from 
human, mouse, camel and bat within 60 minutes. SARS PLpro is also able to fully process 
proISG15 from sheep and shrew within that same time but it has little to no activity against 
jackfish proISG15. In contrast, MERS PLpro is fully capable of processing jackfish proISG15 
but it processes proISG15 from shrew and sheep poorly. MHV PLP2 on the other hand displays 
a much narrower specificity for ISG15s from different species. It shows little to no cleavage of 
proISG15s derived from human, sheep, camel, or bat sources but does show modest activity for 
the northern tree shrew. MHV PLP2 shows strong activity towards jackfish and mouse 
proISG15, the latter activity being consistent with MHV’s pathogenicity as a murine CoV. 

Species specific affinity of ISG15 for SARS-CoV and MERS-CoV PLpros 

Analysis of the cleavage patterns in Figure 2 suggests that SARS and MERS PLpros are 
more promiscuous in recognizing and cleaving various proISG15s compared to MHV PLP2, 
which may relate to the fact that these human pathogens are capable of replicating in different 
hosts, e.g. bats, camels and the shrew, in contrast to MHV which can only replicate in the mouse. 
Differences in cleaving ability of proISG15s by SARS and MERS PLpro may be due to 
sequence and structural differences, which may affect binding affinities for the ISG15s. This 
possibility was explored using Isothermal Titration Calorimetry (ITC) to measure the binding 
affinities of different ISG15s missing the P’-sites (Figure 3). As delSGylases natural substrates 
is an isopetide bond, removal of the P’-sites also allows the avoidance of potential confounding 
factors related to CoV PLPs other prominent function of processing the viral polyprotein to 
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promote viral replication. This part of the process requires that PLPs recognize and cleave the 


peptide bond after the LXGG sequences within the polyprotein, and the ability of PLPs to cleave 
this bond has been shown to be reliant on residues flanking the peptide cleavage site 6 ' 9 . This 
requisite function is not associated with other viral delSGylating proteases such as nairovirus 
vOTUs“ . Therefore, ITC measurements were performed using mature ISG15. Although the 
mature ISG15s resemble the product by the absence of an isopetide bond, binding affinity could 
be assessed independent of any P’-sites contribution associated with using the ISG15-pro-form 
protein. This would also provide a more quantitative understanding of the thermodynamic 
parameters related solely to the species and viral related variances of the protease and ISG15s 
respectively. 

Intriguingly, both SARS-CoV and MERS-CoV PLpro have stronger affinity toward 
hISG15 than shISG15 and mISG15. The dissociation constants (K d ) of SARS-CoV and MERS- 
CoV PLpro for hISG15 are 20.5 ± 4.5 pM and 59.3 ± 4.5 pM, respectively. Unlike their 
affinities for hISG15, the affinities of both PLpros for shISG15 and mISG15 were significantly 
weaker and as a result, a competitive ITC binding assay (Table 2; Figure 3) had to be used to 
determine their K d values. In regards to SARS-CoV PLpro, ITC analyses revealed similar K d ’s 
of 200 ± 41 pM and 198 ± 64 pM for shISG15 and mISG15. More divergence in affinity for 
mISG15 and shISG15 were seen with MERS-CoV PLpro. Like SARS-CoV PLpro, affinity of 
MERS-CoV PLpro for shISG15 was similar with a K d of 147 ± 36 pM. However, a substantially 
weaker affinity of MERS-CoV PLpro for mISG15 was observed (K d of 376 ± 53 pM). 

The measured thermodynamic parameters revealed that both SARS-CoV and MERS- 
CoV PLpro follow a similar trend and prefer hISG15 over shISG15 and mISG15. Specifically, 
the association of shISG15 and mISG15 to PLpros is unfavorable with increasingly higher 
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entropic factors. Although PLpros originating from MERS-CoV and SARS-CoV exhibited a 


similar preference for hISG15, the thermodynamics driving the affinity differ. The SARS-CoV 
and hISG15 binding event was driven by enthalpic factors and was slightly entropically 
unfavorable. In contrast, the MERS-CoV affinity of hISG15 was more balanced possessing an 
entropic and enthalpic component. 

Although some DUB proteases that possess delSGylating activity have been observed to 
engage only the C-terminal domain of ISG15, SARS-CoV has been suggested to interact with 
both domains ’ ' ’ ~ . To gain insight into the relative contributions of the C-terminal domain 
to their full-length counterparts, ITC was performed on utilizing both mouse and human C- 
terminal variants of ISG15. The Kd of 57.6 ± 3.21 pM was observed for ChSIG15. Although the 
binding event still was enthalpically favorable, this contribution was four-fold less that that 
observed for hISG15. Conversely, it possessed a markedly favorable entropic component over 
the slightly disfavorable one of its full-length counterpart. Not unsurprisingly, the trend of 
SARS-CoV having a stronger affinity toward hISG15 over mISG15 was also observed in the C- 
terminal domain of each substrate. The Kd for CmISG15 was considerably weak at 1,870 ± 173 
pM showing marked decreases in both entropic and enthalphic binding contributions suggesting 
that similar to hISG15, mISG15’s N-terminal domain may also be involved in the protease- 
ISG15 binding interface. 

Crystallization of SARS-CoV PLpro bound to the C-terminal domains of human and 
mouse ISG15. 

While X-ray crystal structures have been solved of SARS and MERS PLpro bound to Ub, 
our understanding of the interactions of these enzymes with ISG15 has been largely limited to 
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enzymatic, mutational and computational modeling studies n ’ 15 . To gain molecular insight into 
the specific interactions between PLPs and ISG15, we crystallized and determined the X-ray 
structures of SARS PLpro bound to the C-terminal domains of ISG15 from both human 
(ChISG15) and mouse (CmISG15). Attempts at co-crystallization of full-length ISG15s from 
human and other species with SARS CoV PLpro, as either various covalent adducts or non- 
covalently bound species, were made but were unsuccessful. However, the complex of SARS- 
CoV PLpro bound to the C-terminal domain of human ISG15 (ChISG15), which is the principle 
binding domain of ISG15, crystallized readily 11 . Specifically, we utilized a form of ChISG15 
and CmISG15 that were modified with propargylamine at the c-terminus (ChISG15-PA or 
CmISG15-PA) to form a suicide substrate that covalently links to the active site cysteine. From 
here on in, these covalently modified complexes of SARS PLpro are designated to as SARS- 
CoV PLpro-ChISG15 and SARS-CoV PLpro-CmISG15 for simplicity. 

X-ray SARS-CoV PLpro bound to the C-terminal domain of human ISG15. 

The X-ray structure of the SARS-CoV PLpro-ChISG15 complex was determined to 2.62 

o 

A with two complete copies of SARS-CoV PLpro-ChISG15 within the asymmetric unit (Table 
3). The structure of CoV PLPs contains the classic tertiary fold associated with PLPs consisting 
of the finger, palm, thumb, and ubiquitin-like (Ubl) domains (Figure 4a). However, despite 
sharing the same tertiary fold, the manner in which SARS-CoV PLpro engages the C-terminal 
domain of hISG15 is distinct from that of Ub. Overlaying the SARS-CoV PLpro-ChISG15 
structure with the structure of SARS-CoV PLpro with mono-Ub (PDB entry 4MM3) reveals that, 
compared to Ub, ChISG15 is shifted by approximately 12° in its global orientation (Figure 4b). 
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Closer examination of the binding interface uncovered key differences in how each 


substrate engages the surface of the SARS-CoV PLpro. Although there is only a minor variation 
surrounding a 180° flip of Trpl23 between the two copies of ChISG15 in the asymmetric unit, 
both ChISG15 chains are analogous in how they are accommodated by SARS-CoV PLpro 
compared to Ub (Figure 4e). The differences appear to be driven by how hISG15 and Ub differ 
in their ability to interact with a hydrophobic region consisting of residues Argl67, Met209, and 
Pro248. For Ub, the interactions with this region occur by means of a hydrophobic patch 
consisting of Ile44, Val70, and Leu8 15 . In contrast, ISG15 lacks such a hydrophobic patch. The 
interaction is primarily driven by two distinct sets of hydrophobic interactions: hISG15’s Trpl23 
with Met209 of the protease’s palm region, and hISG15’s Phel49 with a second hydrophobic 
site driven by Pro224 of the protease’s fingers region (Figures 4c, d). This predominate 
hydrophobic binding nature was also seen in the ITC results of SARS-CoV with ChSIG15, with 
an observed increase in entropic favorable interactions when compared with the hISG15 (Table 
2). Surprisingly, while the Ub-bound structure contains other polar and water-mediated 
interactions between Ub and the palm region of the PLpro, such interactions seem to be less 
pronounced in the more charged ChISG15 structure 15 . Overall, there are only a few water 
molecules present within the interface between PLpro and ChISG15 compared to the larger 
number that are observed in the Ub-bound structure potentially indicating a lesser role of water 
mediated hydrogen bond interactions in the binding of hISG15. 

To gain a better understanding of the residues mediating specificity for Ub or ISG15, 
mutants targeting the palm and fingers regions of the of SARS-CoV PLpro were constructed and 
their catalytic activities towards various substrates were determined (Figure 5, Table 1). These 
substrates included Ub- and ISG15-AMC, which are composed of the respective substrate 
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derivatized with the fluorescent leaving group 7-amido-4-methylcoumarin. Additionally, the 


small peptide Z-RLRGG representing the consensus recognition sequence for both Ub and 
ISG15 also attached to AMC was included to probe for potential changes affecting protein- 
protein interactions versus changes that impact intrinsic catalytic activity. 

One of the mutants, M209A, reduces the activity towards Ub-AMC by almost two-fold 
increases its activity towards Z-RLRGG-AMC by nearly five-fold, while retaining wildtype 
levels of activity for ISG15-AMC. From a structural prospective, the reduction of Ub-AMC 
activity is not surprising as M209A shrinks the hydrophobic patch that engages the Ile44 
centered hydrophobic patch in Ub. The observed increase of the activity towards Z-RLRGG- 
AMC is less straightforward. M209 lacks direct interaction with the last five amino acids of Ub, 
or ISG15. However, its replacement by alanine could sterically open up a space that is more 
accommodating for the artificial Z-adduct of the peptide. The other two mutants were found to 
have increased catalytic efficiencies for processing either Ub or ISG15 with a corresponding 
reduction in activity for the other. R167E is over eight times more efficient than the wild type 
enzyme at hydrolyzing Ub but is about 20-times less efficient at hydrolyzing ISG15 (Figure 5c). 
In contrast, Q233E is nearly three-fold more efficient than wildtype at hydrolyzing ISG15 and 
two-fold less efficient hydrolyzing Ub. In the case of Ub, the mono-Ub structure suggests the 
charge flip with R167E may introduce an additional electrostatic interaction with either Arg42 or 
Gln49 of Ub. For the ChISG15 structure, this change occurs in close proximity to the interaction 
between Trpl23 and Met209. Examination of the structure reveals Argl67 may contribute to n- 
n interactions involving Trpl23 and Argl53 of ChISG15 (Figure 4e, Figure 5a). Replacing the 
longer and positively charged arginine with the shorter glutamate removes this contribution. 
Additionally, this charge flip may disrupt the electrostatic interaction between Argl53 and 
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Glul68 of the PLpro, potentially resulting in a loss of affinity for ChISG15. For Q233E, the 


change from a neutral polar to a charged group may create a potential electrostatic repulsion with 
the backbone carbonyl of Ala46 in Ub (Figure 5b). For the ChISG15 structure, in comparison, 
there is not a clear direct interaction between PLpro and ISG15 that would be affected. This 
suggests that the impact of this mutation may stem from internal changes within the PLpro 
regarding the flexibility of the finger region rather than direct interactions with the substrate. 

X-ray Structure of SARS-CoV PLpro bound to the C-terminal domain of mouse ISG15. 

In light of the differences in the nature of interactions between different species’ ISG15s 
revealed by ITC, as well as the apparent plasticity that can occur between Ub and hISG15 within 
the SARS-CoV PLpro active site, insights into the structural sources of this phenomenon were 
sought. To this end, a structure of SARS-CoV PLpro with the C-terminal domain of mISG15 

o 

(CmISG15) was determined to a resolution of 2.4 A (Ligure 6a). The structure of the catalytic 
core domain of SARS-CoV PLpro is consistent with the structure of the ChISG15 SARS-CoV 
PLpro structure. However, a major difference in the orientation of the Ubl domain of SARS-CoV 
PLpro is observed when bound to CmISG15. Unlike previous SARS-CoV PLpro structures 
where the Ubl domains differs only slightly in its position/orientation to the catalytic domain or 
is unobservable because of weak electron density suggesting potential flexibility, the Ubl in the 
SARS-CoV PLpro-CmISG15 structure takes a sharp, almost 90° turn in orientation when 
compared with previous X-ray structures (Ligure 6c). The new orientation does not appear to 
form contacts with the bound CmISG15 suggesting that there is no direct influence on the Ubl to 
adopt this conformation. Lurther analysis of the new orientation reveals a seam made up of 
several electrostatic interactions between the Ubl domain and the thumb domains of the PLpro, 
indicating that this orientation can be stabilized (Ligure 6b). 
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The Ubl domain was not the only global difference between the SARS-CoV PLpro- 


CmISG15 and SARS-CoV PLpro-ChISG15 complexes. Surprisingly, the global orientation of 
CmISG15 is tilted 27° away from the fingers in SARS-CoV PLpro (Figure 6d). Interestingly, 
four amino acids appear to explain the lack of CmISG15’s accommodation in the same 
orientation to that of bound ChISG15. At position 149 in ChISG15 and 147 in CmISG15, there 
is an amino acid difference of a phenylalanine and an isoleucine respectively (Figure 6e). This 
change in CmISG15 eliminates the hydrophobic interaction with the fingers region of SARS- 
CoV PLpro seen in the ChISG15 bound structure and likely aids in the tilt toward the Ubl 
domain. In addition, there is a change from Asn89, Thrl25, and Asnl51 in ChISG15 to Glu87, 
Serl23, and His 149 in CmISG15 (Figure 6f). These amino acid differences create a hydrogen 
bond network between Glu87, Hisl49, and Serl23 in the CmISG15, which is absent in the 
ChISG15 bound structure. This network locks the Glu87 into an unfavorable position for 
CmISG15 to bind in the same manner as ChISG15. The unfavorable position ultimately results 
in repulsion of CmIG15 against PLpro. 

In concert with a loss of the hydrophobic interaction with the finger region and potential 
electrostatic repulsions CmISG15 would incur in binding in the same orientation of ChISG15, 
the CmISG15 bound orientation facilitates the formation of numerous favorable interactions. 
Consistent with the X-ray structures of other delSGylases bound to the C-terminal domain of 
human ISG15 ’ “ , the conserved Trpl21 of mISG15 is centric to the interface. However, 
unlike the interaction observed in the SARS-CoV PLpro-ChISG15 structure, Trpl21 in the 
mISG15 does not insert into a hydrophobic pocket. Instead, mISG15’s Trpl21 and Prol28 form 
significantly smaller hydrophobic interaction with SARS-CoV PLpro via the proteses’s extended 
alkyl chain of Glul68 generated by the hydrogen bond formed between it and Thrl71. In 
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contrast to the bound ChISG15 and outside the last five C-terminus consensus amino acids of Ub 


and ISG15, this weak hydrophobic interaction between CmISG15 and the protease is the only 
one observed. In addition to the hydrophobic binding contributions of mISG15 Trpl21, a water 
molecule is observed to mediate a hydrogen-bonding network between Trpl21’s indole nitrogen 
and several nearby PLpro residues (Figure 6g/h). Including this hydrogen bond network, 
CmISG15 forms almost 40% more hydrogen bonds with the protease then observed with 
ChISG15. One set of these additional electrostatic interactions are centered around SARS-CoV 
PLpro’s Arg83 and CmISG15’s Lysl32. Others are also observed to form between SARS-CoV 
PLpro’s Asnl57, Glnl75 and Argl67 and CmISG15’s Glul30, the carbonyl of Glyl26, and 
His 149 respectively (Figure 6g/h). 

Taken overall, the interface of bound CmISG15 of 872.5 A is similar in size to that of 
the ChISG15 -767 A" interface, but consists of predominately electrostatic interactions whereas 
the former has a significantly greater hydrophobic component. In line with this, a greater 
number of water molecules, approximately 10, can be observed within the interface further 
indicating the degree to which binding may be driven more hydrophilic interactions. This 
characterization of the binding interaction is consistent with the thermodynamic properties 
observed for the binding of CmISG15 to SARS-CoV PLpro (Table 2). Also, when viewed in 
conjunction with the affinity data of mouse and human ISG15s as well as their C-terminal 
domains counterparts suggests that the type of interface formed between CmIS15 and SARS Co- 
V PLpro is less stable then its ChISG15 counterpart and may be one of the contributing factors to 
the overall reduction in binding affinity observed for full-length mISG15 when compared to 
hISG15 via ITC. 
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X-ray structure of the full-length, unbound form of mouse ISG15 and comparison human 
ISG15 

Whereas some viral delSGylating enzymes, such as vOTUs, are thought to exclusively 
interact with the C-terminus of ISG15, SARS PLpro has been shown to contain two Ub-binding 
sites that likely engage the N-terminus of ISG15s u ’ 15 . With only the full-length hISG15 
structure available, the impact of ISG15 biodiversity within the N-terminal region has been 
difficult to assess. To address this issue, the X-ray structure of mISG15 with both domains was 
determined to assess conformational variability of mISG15 compared to hISG15. Initially, only 

o 

a low resolution (4 A) structure could be determined. Truncation of the last five non-structured 

o 

amino acids led to a structure that could be determined to a higher resolution of 3.25 A. Using 
the program Define Secondary Structure of Proteins (DSSP), the secondary structure of mISG15 
was determined to be comparable to hISG15. mISG15 contains nine beta sheets, two helices, but 
only two 3 10 helices within the C-terminal domain; this differs from the hISG15 which contains 
two 3 io helices per domain of ISG15 (Figure 7a) . The asymmetric unit consists of 2 twisting 
filaments, containing 10 copies of mISG15. Upon examining the differences between the 
monomers within the asymmetric unit, some flexibility between N- and C-terminal domains was 
observed (Figure 7b/c). This flexibility was probed by aligning the C-terminal domains of each 
mISG15 monomer and measuring the angle between the point of divergence and the farthest N- 
termini. There is a 20.2° range of motion between the N- and C-terminal domains of mISG15 
(Figure 7b). Comparable flexibility was also observed in hISG15 structures, with a 18.3° range 
of motion (Figure 7c) (PDB: 1Z2M, 3R66, 3PSE) 27;29;30 . 

While each individual domain of hISG15 and mISG15 are similar in secondary structure 
and both show some structural variability between monomers, there is a drastic difference in the 


17 




overall conformations of the tertiary structures between mISG15 and that of hISG15. This 


surprising difference is observed in the overall arrangement of the N- and C-terminal domains. 
Specifically, the twist about the C- and N-terminal domain of mISG15 in relation to hISG15 
ranges from 43.0° to 66.9° (Figure 8a). Closer investigation reveals that differences in tertiary 
arrangement can be attributed to several molecular interactions within the structures driven by 
the primary sequence differences of mISG15 compared to hISG15. Specifically, the presence of 
Asp79 creates a kink in the hinge region of hISG15 as a result of the carboxylate group of Asp79 
forming a hydrogen bond with the hydroxyl group of ThrlOl (Figure 8b). This interaction does 
not occur in mISG15 since Ser77 replaces Asp79. 

While amino acids in the hinge may play a part in the different conformational trajectory 
allowed in the region between the two domains in mISG15, the twisted feature of mISG15 may 
be stabilized by amino acid interactions in the core region of the protein. This core region is 
centered near Phe41, at the interface between the two domains. The hydrophobic interaction 
between the N-terminal and C-terminal domain at Leul34 and Phe41 may stabilize mISG15s’ 
twisted arrangement (Figure 8b). Other residues in close proximity, such as Pro39, appear to 
further stabilize this hydrophobic interaction in mISG15. Specifically, Pro39 in mISG15 forms 
hydrophobic interactions with Phe41 contributing to the stabilization of mISG15s’ tertiary 
arrangement (Figure 8b). 

In addition to the presence of favorable interactions, steric clashes are also likely 
responsible for preventing the occurrence of a shared tertiary arrangement between the mISG15 
and hISG15 species. When the N-terminal and C-terminal domains of mISG15 are configured to 
the arrangement of hISG15, there are no obvious clashes that prohibit mISG15 from adopting 
this conformation (Figure 8c). However, hISG15 may be incapable of adopting the mouse 
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conformation due to the potential steric clash between His39 and Glul39 that typically form a 


water mediated interaction (Figure 8c). Thus, the presence of the His39 residue may prevent 
hISG15 from configuring to the mISG15 conformation. 

Beyond the effects on tertiary structure, the sequence divergence between mISG15 and 
hISG15 also impacts the potential binding surfaces. These differences in amino acid sequence 
also give rise to different electrostatic potential plots. Unlike mISG15, hISG15 contains an 
uninterrupted negative surface band spanning across the two domains (Figure 8d). When the 
surface of mISG15 corresponding to the same area is examined, this negative band dissipates as 
the comparable area on mISG15 is scattered with positive and neutral areas. 

DISCUSSION 

Differences in hISG15 and mISG15 structure influence recognition by PLPs 

Among Ub and Ub-like proteins, ISG15 is unique for more than its divergent amino acid 
sequence. Similar to other Ub-like proteins, ISG15 has not been observed to form polymeric 

T 1 

chains like Ub . However, ISG15 is the sole family member of ubiquitin-like modifiers that is 
comprised of two Ub-like folded domains. These two domains are tethered by a polypeptide 
linker that would suggest that the two domains can move freely and independent of each other. 
However, the four structures of hISG15 alone, or bound to viral proteins suggest the opposite 
(Figure 7c) ’ ’ . Akin to polymeric Ub, ISG15 appears to have relatively limited 

conformations between the two domains. The structure of mISG15 (Figure 7b) furthers this 
assertion with one important caveat, one species ISG15s domain configuration may not be 
necessarily representative of others, or ISG15s in general. 
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In both hISG15 and mISG15 there is a hydrophobic interface mediated by a conserved 


phenylalanine (Phe41, see Figure 8). The highly-conserved nature of this phenylalanine and 
other surrounding residues suggests that ISG15s in general likely utili z e hydrophobic forces to 
adhere the two domains together (Figure lb, residues involved are boxed in blue). However, as 
shown with mISG15, the influence of the core on orientation of the domains can vary. In 
mISG15, residues forming the core solely dictate its inter-domain orientation, whereas in hISG15 
additional electrostatic interactions involving Glul39 and His39, and to a lesser extent ThrlOl 
and Asp79, appear to play an additional role to favor one orientation over another (Figure 8b,c). 
Intriguingly, the pairing of Glul39 and His39 is extremely unique to hISG15. Although Glul39 
or another acidic residue is well conserved at that position, His39 is typically a proline. In 
hISG15, it appears that this favorable electrostatic interaction may promote hISG15’s domain 
orientations. Multiple crystal structures reveal this conformation to be consistent despite 
differences in space groups and crystallization conditions, suggesting the observed tertiary 
structure a stable and likely preferred conformation (Figure 7c). In addition to favoring 
hISG15’s conformation, His39 may also act as a steric block preventing hISG15 from adopting a 
similar tertiary structure and conformation as that of mISG15. 

The observation of a potential steric block may result in MHV PLP2 not being able to 
recognize hISG15 which is why no cleavage of pro-hISG15 is observed (Figure 2). In contrast, 
no apparent steric hurdle is present for mISG15 in adopting a hISG15 like inter-domain 
orientation which maybe why SARS and MERS PLpro are able to readily recognize mISG15. 
The lack of a steric hurdle might suggest that there are more allowable domain arrangements of 
mISG15 than hISG15 in solution that allow for a broader spectrum of PLPs to recognize 
mISG15. The potential uniqueness of h!SG15’s structure may also fall in line with the inability 
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of other viral proteins, such as influenza NS1, to effectively engage ISG 15s beyond those of 

on 

human and primates . Moreover, the distinctiveness of hISG15 may also fit into the recent 
assertions that hISG15 plays a divergent immune regulation role in humans than in other animals 

32 

Viral USP-Like PLPs accommodation of Ub versus ISG15s 

Ub interactions with proteins possessing an ubiquitin binding site have been observed to 
be reliant largely on the involvement of a hydrophobic patch comprised of several residues 
surrounding Ile44 15 . The absence of such a comparable patch in ISG15 as well as the more 
generally charged surface of hISG15 naturally spurred speculation on mechanisms behind how 
PLPs may engage ISG15s. The SARS-CoV PLpro-ChISG15 complex reveals that in PLPs the 
binding interface may have evolved to recognize specific features of ISG15s from different 
species outside the five-amino acid, C-terminal sequence (Arg-Leu-Arg-Gly-Gly), they share 
with Ub. For hISG15, this includes not only the highly conserved Trpl23 of ISG15s, but also 
specific interactions that are unique to hISG15. For instance, Arg87 in hISG15 is often 
substituted for a lysine residue that is too short form any water mediated interactions with the 
SARS-CoV PLpro’s Asp230 and Ser222. Interestingly, Phel49, which forms a second 
hydrophobic interaction site with the protease’s Pro224, is one of the three amino acids recently 
implicated in species specificity among nairovirus vOTUs . Residue changes between hISG15 
and mISG15 at the other two positions, Asnl51 and Asn89, appear to impede the 
accommodation of mISG15 in the active site of SARS-CoV PLpro compared to the more 
favorable hISG15 orientation. This appears to advocate that differences at these three ISG15 
positions have a broader range of influence beyond only one class of viral proteases and could 
conceivable represent an evolutionary pressure that underlies part of ISG15 sequence diversity. 
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Also, SARS-CoV residues that interact with these residues as well as SARS-CoV PLpro 


deubiquitinating and delSGylating altering mutants, R167E, M209A, and Q233E, may offer 
advantageous starting points for developing SARS-CoV PLpros with directed shifts in substrate 
specificities. 

Beyond the C-terminal domains interaction with PLPs, the N-terminus of human ISG15 
has been previous proposed to interact with a ridge helix spanning the conical PLP thumb 
domain with the Ubl domain to enhance affinity. Currently, no X-ray structures of an PLP with 
full-length ISG15 exist and previous computational models utilizing existing SARS-CoV PLpro 
structures bound with mono-Ub have had difficulty reconciling the significant gap that occurs 
between hISG15’s N-terminal domain and the protease when anchoring C-terminal hISG15 
domain on bound Ub. The SARS-CoV PLpro-CxISG15 structures offer two synergistic 

o 

explanations. First, the 6.3 A shift of ChISG15 relative to bound Ub translates the N-terminal 
domain towards the ridge helix (Figure 9a). In addition, the alternate Ubl conformation found in 
the SARS-CoV PLpro-CmISG15 structure reveals that such a conformation translates the ridge 

o 

helix up to 14 A toward the location of the ridge helix (Figure 9b/c). Interestingly, comparing 
this model to the X-ray structure of SARS-CoV PLpro bound to K48-linked di-Ub highlights that 
different facets of the Ub-like fold found in the N-terminal domain ISG15 are likely involved in 
the interaction (Figure 9d). Specifically, in this model a triple serine repeat, Ser20-22, and 
Glu27 located in hISG15’s P2a3 loop and a3 respectively point toward the protease’s ridge 
helix creating a surface that is available to engaged by the bevy of charged and polar residues 
located on SAR-CoV PLpro’s ridge helix that has been previously implicated in binding 15 . 
These potential electrostatic interactions may contribute to the four times larger enthalpy binding 
component observed in the interaction of full-length over the C-terminal domain of hISG15 
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when being accommodated by SARS-CoV PLpro. The impact of mISG15’s divergent domain 


orientation from hISG15 is also apparent. Initial molecular modeling of the hISG15 structure 
onto the CmISG15-bound SARS-CoV PLpro structure with CmISG15 as an anchor reveals a 
steric clash with the ridge helix. This is the same for hISG15 drawn from any of its known X-ray 
structures. However, the divergent inter-domain orientation found in full-length mISG15 
structure determined here allows for the mISG15’s N-terminal domain to fit unobstructed (Figure 
9e). The significant domain-domain orientations differences between hISG15 and mISG15 
results in a different facet of mISG15 facing the protease. Specifically, mISG15’s al[13 loop 
comprises the surface oriented towards the ridge helix presenting a polar interface that ITC 
suggests may be involved in forming additional electrostatic interactions beyond those of the C- 
terminal domain. However, when comparing the ITC data for full-length and its CmISG15 
counterpart, these additional interactions might potentially come at an ethalphic cost. This may 
suggest that different ISG15 N-terminal domains may engage the protease’s ridge helix to 
differing degrees and thermodynamic characteristics. The need for PLPs not only to 
accommodate surface differences between species ISG15s, but also divergent inter-domain 
orientations highlight the difficulties for a CoV PLP to be active to all species ISG15s. This 
could also be perceived as a possible benefit to ISG15’s unique tandem Ub-like arrangement. 
Also, the improved structural perspective on how the N-terminus of ISG15s from different 
species may engage CoV PLPs provides a clearer path towards utilizing this region to influence 
delGylating activities of CoV PLPs. As a result, combining this information along with the 
alterations possible in the C-terminal domain, fresh tools to addressing the role of 
deubiqutination and delSGylation through the use of reverse genetics systems can be envisioned. 

Possible evolution of CoV PLP recognition of species variances in ISG15s 
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Overall, the X-ray structural, enzymatic and biophysical data point to CoV PLP 


delSGylase activities being sensitive to species specific amino acid differences within ISG15. 
Intriguingly, SARS-CoV and MERS-CoV PLpros, whose viruses replicate in a wide range of 
hosts, recognize and cleave proISG15 from almost all of the species tested. In contrast, the 
mouse specific MHV PLP2 is limited predominantly to the mouse substrate. This wider range of 
specificities appears substantially larger than that found recently in nairovirus vOTUs. This may 
imply that with the greater binding interface provided by the palm, fingers and thumb domains of 
the ubiquitin specific protease fold, CoV PLPs can either engage a wider array of ISG 15s then 
that of vOTUs, or be highly selective for just one or two ISG15s as in the case of MHV. 

Intriguingly, in all of the PLPs examined, their potential ability to engage different 
species ISG15s is not fully restricted to only those from species their parent viruses infect. 
Understandably, this could be a by-product in the evolutionary process of a specific viral 
protease seeking to optimize towards a certain species ISG15 and inadvertently picking up 
enzymatic activity towards another ISG15. Or, in an environment where a virus only has to 
optimize its replication in one species, there is less selective pressure and hence a greater chance 
of losing the ability of recognizing ISG15s from other species. Alternatively, some of these types 
of off-species PLP activities could be indicative of evolutionary memory for ISG15s. In essence, 
this could give a possible view into the zoonotic evolutionary history, or potential future 
zoonotic drift for a certain virus. Naturally, a wider sampling of CoV PLPs affinities for certain 
species ISG15s, knowledge of what species their parent viruses infect, and a multitude of reverse 
genetics experiments will be necessary to discern which of the above scenarios in taking place in 
the evolution of virus recognition of host ISG15s. 

MATERIALS AND METHODS 
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Construct, Expression and Purification of PLPs for the ISG15 Protease Activity Assay and 
ITC. 

MHV PLP2 was expressed and purified as previously described ~ ’ . SARS-CoV PLpro 
in expression vector pET21a and MERS-CoV PLpro in pET15b were transformed into 
Escherichia coli BL21(DE3) competent cells (New England Biolabs) by heat shock. Cells were 
grown at 37°C in LB broth supplemented with 100 pg/mL of ampicillin to ODeoo 0.6-0.8 and 
expression induced with 0.8 mM isopropyl-[3-D-thiogalactopyranoside (IPTG) at 25°C overnight 
for SARS-CoV PLpro and 1 mM IPTG at 18°C overnight for MERS-CoV PLpro. Cells were 
collected by centrifugation at 6000xg for 10 minutes and stored at -80°C. Cells were lysed in 
Buffer A (20 mM Tris [pH 7.5 for SARS-CoV PLpro, pH 7.0 for MERS-CoV PLpro], 500 mM 
NaCl, 10 mM |3-mercaptoethanol [BME]) supplemented with lysozyme for 30 minutes at 4°C, 
followed by sonication on ice at 50% power with a 50% duty cycle for a total of 6 minutes. 
Insoluble protein was removed by centrifugation at 70,600xg for 30 minutes and the supernatant 
filtered through a 0.80 pm filter. The clarified supernatant was flowed over high density nickel 
agarose beads (Gold Biotechnology, Olivette, MO) pre-equilibrated with Buffer A. The column 
was washed with 10 column volumes of Buffer A supplemented with 30 mM imidazole, and the 
protein eluted with 10 column volumes of Buffer A supplemented with 300 mM imidazole. The 
PLpro was further purified by size exclusion chromatography using a Superdex 200 column (GE 
Healthcare, Pittsburgh, PA) pre-equilibrated with 100 mM NaCl, 5 mM HEPES [pH 7.5 for 
SARS-CoV PLpro, pH 7.0 for MERS-CoV PLpro] and 2 mM DTT. 

Purification of SARS-CoV PLpro for Complexation and Crystallization. 

Purification of SARS-CoV PLpro in expression vector pETlla for complexing with the 
C-terminal domain of human ISG15 (ChISG15) and mouse ISG15 (CmISG15) was adapted from 
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the previously described method . The cells were chemically lysed by resuspending them in 
150 mL of Buffer B (20 mM Tris [pH 7.5], 10 mM BME) and lysozyme and incubated at 4°C for 
30 minutes. The suspension was sonicated on ice at 50% power with a 50% duty cycle for a total 
of 6 minutes and centrifuged for 30 minutes at 40,900xg. The cell lysate was filtered and 
subjected to a 40% ammonium sulfate fractionation then centrifuged again for 30 minutes at 
40,900xg. The resulting pellet was resuspended in 250 mL of 1 M ammonium sulfate, 20 mM 
Tris [pH 7.5], and 10 mM BME and incubated at room temperature for 1 hour. The suspension 
was filtered and loaded onto a 50 mL Phenyl-Sepharose CL-4B column (GE Healthcare, 
Pittsburgh, PA) equilibrated with 1.5 M ammonium sulfate, 20 mM Tris [pH 7.5], and 10 mM 
BME. The protein was eluted using a 10-column-volume gradient to 100% Buffer B and washed 
with 2 additional column-volumes of 100% Buffer B. The fractions were pooled together and 
then diluted fivefold with Buffer B. The protein was loaded onto a MonoQ 10/100 column (GE 
Healthcare, Pittsburgh, PA) equilibrated with Buffer B. The protein was eluted using a 10- 
column-volume gradient to 100% of a buffer composed of 0.5M NaCl, 20 mM Tris [pH 7.5], and 
10 mM BME with the initial flow-through off the column collected. The flow through was 
concentrated and put into dialysis in a 50 mM NaCl, 20 mM Tris [pH 7.5], and 10 mM BME 
buffer at 4°C overnight. 

Construct, Expression and Purification of proISG15s and mature ISG15s. 

ISG15s from human ( Homo sapiens ; Accession: AAH09507.1), mouse (Mus musculus ; 
Accession: AAB02697.1), northern tree shrew ( Tupaia belangeri ; Accession: AFH66859.1), 
sheep (Ovis aries; Accession: AF152103.1), dromedary camel ( Camelus dromedarius; 
Accession: XP_010997700.1), vesper bat ( Myotis davidii ; Accession: ELK23605.1), and 
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jackknife fish (Oplegnathus fasciatus\ Accession: BAJ16365.1) in both pro and mature forms 
were prepared as described elsewhere 23 . 


ISG15 Protease Activity Assay. 

Activity assays of SARS-CoV PLpro, MERS-CoV PLpro, and MHV PLP2 with purified 
northern tree shrew proISG15 (pro-nsISG15), sheep proISG15 (pro-shISG15), fish proISG15 
(pro-fISG15), mouse proISG15 (pro-mISG15), camel proISG15 (pro-cISG15), bat proISG15 
(pro-bISG15), and human proISG15 (pro-hISG15) were adapted from the previously described 
methods 23 . 

ITC of ISG15 with PLpros from MERS-CoV and SARS-CoV. 

ITC was performed using a Microcal PEAQ-ITC (Malvern, Worcestershire, UK). There 
were 19 injections of 2 pL each at 25°C with a reference power of 6 pcal/s. The mature forms of 
ISG15s along with PLpros from MERS-CoV and SARS-CoV were dialyzed at 4°C in 50 mM 
HEPES [pH 7.4], 200 mM NaCl, and 1 mM DTT. All experiments were run in duplicate. For 
direct binding experiments, 227 pM and 276 pM of SARS-CoV and MERS-CoV PLpro 
respectively were placed in the cell with 2.3-2.6 mM of mature hISG15 in the syringe. For direct 
binding experiments of SARS-CoV with ChISG15, 303 pM was placed in the cell with 3.4 mM 
in the syringe respectively and for SARS-CoV with CmISG15, 393 pM was placed in the cell 
with 3.8 mM in the syringe respectively. For competitive experiments related to SARS-CoV 
PLpro, mixtures containing 100 pM of protease with 50 pM of either sheep ISG15 (shISG15), or 
mISG15 was placed in the cell with 1 mM of mature hISG15 in the syringe. For competitive 
experiments related to MERS-CoV PLpro, mixtures containing 220 pM and 270 pM of protease 
with 110 pM and 170 pM of shISG15 and mISG15 respectively was placed in the cell. The 
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syringe contained 2.3 mM and 2.7 mM of mature hISG15 in the syringe for shISG15 and 


mISG15 assays respectively. The data were processed using Microcal PEAQ-ITC Analysis 
Software. 

Functional studies of SARS-CoV PLpro Mutants. 

The SARS-CoV pET-15b-PLpro mutants (residues 1541-1855 of the SARS-CoV viral 
polyprotein) were generated using site-directed mutagenesis and the QuickChange® approach 
(Agilent). Expression and purification for the wild type and each mutant of the SARS-CoV 
PLpro were performed as previously described l4 . 

The steady-state kinetic parameters of SARS-CoV PLpro wild-type and mutants were 
determined for three different ubiquitin-based fluorescent substrates, utilizing 7-amino-4- 
methylcoumarin (AMC), commonly used to assess the protease, deubiquitinating, and 
delSGylating activity of PLPs, including a small peptide substrate, Z-RLRGG-AMC (Bachem), 
Ub-AMC (LifeSensors, Inc.), and ISG15-AMC (Boston Biochem/R&D Systems). Kinetic 
assays with Ub-AMC and ISG15-AMC were performed on the same day and side-by-side in the 
same assay plate to directly compare the enzymatic activity of SARS-CoV PLpro to that of each 
of the mutants. The steady-state kinetic studies were also repeated for the wild-type and mutants 
approximately five months apart and the resulting duplicate data were combined for analysis. 
Kinetic assays with the peptide substrate were also performed in triplicate. For all experiments, 
the assay conditions, i.e. buffering conditions and assay volume etc, were setup as previously 
described “ . The exception was that the stock substrates purchased from the vendors had 
different lot numbers. The steady-state kinetic data obtained from separate experiments 
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performed on different days and with different substrates lot numbers helped to ensure that the 


trends in the resulting kinetic parameters were reproducible. 

The enzymatic activity of PLpro-mediated hydrolysis of the fluorophore, AMC group 
was determined using a BioTEK Synergy HI multimode microplate reader at 25°C with a 
wavelength of excitation at 360 nm (bandwidth=40 nm) and an emission wavelength of 460 nm 
(bandwidth=40 nm). The change in the relative fluorescence as a function of time (RFU/min) 
was monitored over a sufficient time period to allow determination of the enzymatic rate in the 
steady-state region. For the ISG15-AMC assay, the substrate concentrations were varied from 
0.2 pM up to 19.2 pM. The reactions were initiated by the addition of enzyme with the final 
enzyme concentrations as follows: 0.48 nM WT, 0.23 nM Q233E, 0.23 nM M209A, or 7.3 nM 
R167E. For the Ub-AMC assay, substrate concentrations were varied from 0.5 pM to 17.6 pM. 
The final enzyme concentrations were 3.7 nM WT, 7.3 nM Q233E, 7.3 nM M209A, or 0.23 nM 
R167E. For the Z-RFRGG-AMC assay, the concentrations of substrate were varied from 0.8 pM 
to 50 pM and the final concentration of the wild type enzyme was 0.14 pM. To capture the 
initial rate of peptide hydrolysis for the M209A mutant a lower enzyme concentration of 25 nM 
was used. As is typical for SARS-CoV PFpro, the enzyme could not be saturated with the Ub- 
AMC and Z-RFRGG-AMC substrates. As such, the kinetic response of the enzyme to these 
substrates was linear and thus the data were fit to a line to approximate the catalytic efficiency 
(k cat /K m ) for each enzyme. For the ISG15-AMC assays, the data were fit to the Michaelis- 
Menten equation to determine the associated kinetic parameters (k cat , K m , and k cat /K m ) for each 
enzyme “ . Saturation was not attained with the R167E mutant enzyme for ISG15-AMC and 
therefore this kinetic data was also fit to a line to determine the apparent (k cat /K m ). The errors 
associated with each kinetic parameter were obtained from the best-fit line or curves for each 
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mutant. All data, from separate experiments, were included in the fits to arrive at the final errors 


(Table 1). 

SARS-CoV PLpro-CmISG15 and SARS-CoV PLpro-ChISG15 Complex Formation. 

Expression of CmISG15, or ChISG15 occurred using a vector pTYB2 and was purified 
as previously described to form a propargylamine-derivatized thioester product (CmISG15-PA, 
ChISG15-PA) “ . Briefly, to obtain complex, purified protease was added directly to the mixture 
in equi-molar ratios, and incubated for 2-4 hours at room temperature and left 4°C overnight. To 
further purify the complex, anion exchange chromatography was used, eluting from a MonoQ 
10/100 column using a linear gradient from 0 to 1 M NaCl (SARS-CoV PLpro-CmISG15) or 
250 mM NaCl (SARS-CoV PLpro-ChISG15) with 50 mM Tris (pH 8.0 for SARS-CoV PLpro- 
CmISG15 and pH 9.0 for SARS-CoV PLpro-ChISG15), followed by size exclusion 
chromatography on a Superdex 75 column (GE Healthcare, Pittsburgh, PA) pre-equilibrated with 
100 mM NaCl, 5 mM HEPES [pH 7.5], and 2 mM DTT. For SARS-CoV PLpro-ChISG15, an 
additional purification step prior to anion exchange chromatography was used to eliminate 
residual ChISG15 by size exclusion chromatography using a Superdex 200 column (GE 
Healthcare, Pittsburgh, PA) pre-equilibrated with 100 mM NaCl, 50 mM Tris [pH 7.5], 2 mM 
DTT. 

Crystallization of SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and 
mISG15. 

Purified SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and mISG15 were 
screened against a series of Qiagen NeXtal suites by hanging drop using a TTP Labtech 
Mosquito (TTP Labtech, Herfordshire, UK) at 8.8 mg/ml, 8.88 mg/ml, and 16 mg/ml 
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respectively. For SARS-CoV PLpro-CmISG15, the initial screen yielded the best crystals in a 


solution containing 65% (vol/vol) MPD and 0.1 M Tris (pH 8.0). These crystals were optimized 
using the Additive HT Screen from Hampton Research. The final SARS-CoV PLpro-CmISG15 
crystals were obtained through vapor diffusion using a 500 pL reservoir with 4 pL hanging drops 
mixed 1:1 with protein solution and reservoir solution, which also contained 0.25 pi of 30% 
(w/v) Trimethylamine N-oxide dihydrate. For SARS-CoV PLpro-ChISG15, the initial screen 
yielded the best crystals in a solution containing 0.2 M lithium sulfate, 0.1 M Bis-Tris (pH 6.5) 
and 25% (w/v) PEG3350. The initial crystal conditions for SARS-CoV PLpro-ChISG15 crystals 
were optimized along salt and PEG3350 gradients as well as using the Additive Screen from 
Hampton Research. The final SARS-CoV PLpro-ChISG15 crystals were obtained through vapor 
diffusion using a 500 pL reservoir with 4 pL hanging drops mixed 1:1 with protein solution (6.99 
mg/mL) and reservoir solution (0.1 M lithium sulfate, 0.1 M Bis-Tris, [pH 6.5], 22% PEG3350), 
which also contained 0.25 pi of 30% (v/v) glycerol. For mISG15, the initial screen yielded the 
best crystals in a solution containing 0.2 M ammonium sulfate, 0.1 M tri-sodium citrate pH 5.6, 
15% (w/v) PEG4000. This condition was further optimized along buffer, pH, and PEG4000 
gradients in addition to using the Additive Screen from Hampton Research. In conjunction with 
using these optimization methods, the mISG15 was shorted by five amino acids. The final 
crystals of the five amino acid shortened mISG15 were obtained through vapor diffusion using a 
500 pL reservoir with 4 pL hanging drops mixed 1:1 with protein solution (6.9 mg/mL) and 
reservoir solution (0.2 M ammonium sulfate, 0.1 M sodium acetate pH 4.6, 12% PEG4000 and 
0.2 M sodium malonate). 

Crystals of SARS-CoV PLpro-CmISG15, SARS-CoV PLpro-ChISG15, and mISG15 
were collected and flash frozen in liquid N 2 . Cryogenic solutions for SARS-CoV PLpro- 
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CmISG15, SARS-CoV PLpro-ChISG15 constituted their respective mother liquors. For 


mISG15 crystals, they were passed from a 5% to a 12% solution 1:1:1 of glycerol, dimethyl 
sulphoxide (DMSO), and polyethylene glycol known as EDG 35 . Data sets were collected at the 
Advanced Photon Source (Argonne National Labs, Argonne, IL). A data set for SARS-CoV 
PLpro-CmISG15 was collected at the LS-CAT beamline 21G at a wavelength of 0.9786 A using 
a MAR300 detector, whereas data sets for SARS-CoV PLpro-ChISG15 and mISG15 were 
collected at the SER-CAT beamlines 22ID and 22BM at 1 A using MAR300hs detectors. All 
data sets were collected at 100 K. 

Data Processing and Structure Solutions. 

Data sets were indexed, integrated and scaled using HKL-2000 36 . All the structures were 
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solved by molecular replacement using Phaser . Subsequently, each structure was rebuilt 
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initially using Autobuild followed by successive rounds of manual model building and 
refinement using Coot 39 and Phenix 40 . The initial solution for the SARS-CoV PLpro-CmISG15 
complex was achieved by using the catalytic core of a previous SARS-CoV PLpro structure 
(PDB entry 3E9S). Density from the last 10 amino acids of the CmISG15 molecule served as an 
anchor for the initial global placement of the CmISG15 from a Erve nairovirus vOTU-CmISG15 
complex (PDB entry 5JZE). This partial model was then used as a search model along with the 
Ub-like SARS-CoV PLpro domain from PDB entry 4MM3 to obtain a complete global model 
using Phaser . For the SARS-CoV PLpro-ChISG15 complex, an initial molecular replacement 
solution was obtained by using the core and other elements from the SARS-CoV PLpro- 
CmISG15 as a search model. A partial molecular replacement solution using Phaser for 
mISG15 was obtained by searching with the CmISG15 from the Erve nairovirus vOTU- 
CmISG15 complex structure (PDB entry 5JZE). This partial model was used in a sequential 
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'J'l 

Phaser run using a mISG15 N-terminal domain homology model based on the previously 


solved hISG15 structure (PDB entry 1Z2M). All structures were validated using Molprobity 41 
and have good Ramachandran statistics: SARS-CoV PLpro-CmISG15 (96% favored and 4% 
allowed), SARS-CoV PLpro-ChISG15 (96% favored and 4% allowed) and mISG15 (99.08% 
favored and 0.92% allowed). All structures have been deposited in the protein data bank. Codes 
can be found on Table 3. 

Electropotential plots. 

Figure renderings involving electropotential plots were performed using the PDB2PQR 

42 

server and the surface generated using the adaptive Poisson-Boltzmann solver (APBS) . 

ACCESSION NUMBERS 

All structures have been deposited in the Protein Data Bank (PDB code 5TL7 for SARS-CoV 
PLpro-CmlSG15, 5TL6 for SARS-CoV PLpro-ChISG15, and 5TLA for mISG15). Source sequences for 
ISG15s were GenBank numbers AAH09507.1 for Homo sapiens, AAB02697.1 for Mus musculus, 
AFH66859.1 for Tupaia belangeri, ELK23605.1 for Myotis davidii, BAJ16365.1 for Oplegnathus 
fasciatus, and NCBI Reference Sequence XP 010997700.1 for Camelus dromedarius. Source sequences 
for PLPs were GenBank numbers AFS88944.1 for MERS CoV PLpro, AKP80587.1 for PEDV PLP2, 
AHM88399.1 for PDCoV PLpro, and UniProtKB/Swiss-Prot numbers P0C6U8 for SARS CoV PLpro 
and P0C6V0 for MHV PLP2. 
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Figure and Table Legends 


Table 1. Kinetic Parameters of SARS-CoV PLpro WT and mutants with different ubiquitin- 
based fluorescent substrates. 

Table 2. Isothermal Titration Calorimetry. 

Table 3. Data collection and refinement statistics. 

Figure 1 (a) The PLPs from SARS CoV (Accession number P0C6U8), MERS CoV (Accession 
number AFS88944.1), MHV (Accession number P0C6V0), PEDV (Accession number 
AKP80587.1), and Porcine deltacoronavirus (PDCoV) (Accession number AHM88399.1). The 
sequence numbering ruler is based on the SARS-CoV PLpro sequence. Residues forming the 
catalytic triad are boxed in purple. The secondary structure of SARS-CoV PLpro based on 
DSSP is denoted, with a- and 3io-helices in reddish orange and P-sheets in wheat. Residues 
forming the “ridge” of SARS-CoV PLpro are marked with an asterisk (*). Residues of SARS- 
CoV PLpro mutated in this study are denoted (“x”). (b) Sequence alignment of ISG15s from 
mouse ( Mus musculus ; Accession: AAB02697.1), human ( Homo sapiens ; Accession: 
AAH09507.1), dromedary camel Camelus dromedarius; Accession: XP_010997700.1), vesper 
bat ( Myotis davidii; Accession: ELK23605.1), sheep (Ovis aries; Accession: AF152103.1), 
northern tree shrew ( Tupaia belangeri; Accession: AFH66859.1), and jackknife fish 
(Oplegnathus fasciatus'. Accession: BAJ16365.1) along with ubiquitin. The sequence numbering 
ruler is based on mouse ISG15. Secondary structure of mouse ISG15 based on DSSP is denoted 
in purple. The hinge region between P4 and P5 is denoted using a dashed line in teal. 
Percentages indicate sequence identity relative to mouse ISG15. Residue positions identified in 
this study that form key interactions with SARS-CoV PLpro are boxed in red. Residue positions 
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implicated in driving the specific tertiary fold features of mouse and human ISG15 are boxed in 
blue. 

Figure 2 Cleavage profiles of SARS-CoV PLpro, MERS-CoV PLpro, and MHV PLP2 with 
proISG15s derived from human, sheep, northern tree shrew, jackknife fish, mouse, dromedary 
camel, and vesper bat. Reactions were conducted at 37°C with 10 pM of each proISG15 and 20 
nM of each PLpro/PLP2 with samples taken at the indicated time points. 

Figure 3 ITC binding isotherms with the raw heat (top panel) and the integrated heats of 
injection (bottom panel) shown for each interaction. 

Figure 4 (a) Cartoon representation of the SARS-CoV PLpro complexed with ChISG15 (blue). 
Secondary structure of the PLpro is designated, with helices and loops rendered in reddish 
orange and P-sheets rendered in wheat. The structural domains of the PLpro are identified, 
consisting of the fingers (white), palm (purple), thumb (green), and Ubl (red) labeled with 
colored boxes. The vinyl thioether propargylamine linker is colored green, and the Zn(II) ion in 
light purple, (b) Overlay of ChISG15 (blue) bound to SARS PLpro (reddish orange) compared 
to a Ub (light brown)-bound structure (PDB entry 4MM3) based on a secondary structure 
alignment of the respective PLpros. The approximate degree of shift in orientation of analogous 
a-helices is indicated, (c) Comparison of the interaction of ChISG15 (blue) versus Ub (light 
brown) with the hydrophobic patch of SARS-CoV PLpro, with the site of interaction within the 
overall structure indicated by an arrow. Side chain conformations of the PLpro are colored 
orange for the ChISG15 bound structure, gray for the Ub-bound structure. Residues E127-D133 
of ChISG15 and A46-T55 of Ub were removed for clarity. Inter- and intramolecular distances 

o 

(A) are shown by black dashes, (d) Comparison of ChISG15 versus Ub at the site of the 
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additional hydrophobic interaction in the ChISG15-bound structure. Colored as in (c). 

o 

Intermolecular distances (A) are shown as as black dashes, (e) Interchain variability of ChISG15 
Trpl23 in binding to SARS-CoV PLpro. The two copies within the asymmetric unit were 
overlaid based on a secondary structure alignment of the PLpros. Side chains of the PLpro 
colored in orange correspond to the ChISG15 colored light blue, and the PLpro side chains 
colored burgundy to the ChISG15 colored a darker blue. Intermolecular distances are colored in 
black. Residues E127-D133 of ChISG15 omitted for clarity. 

Figure 5 (a) and (b) Sites within the ChISG15- and Ub-bound SARS-CoV PLpro structures 
corresponding to the mutations causing increased Ub-AMC activity (top) versus increased 
ISG15-AMC activity (bottom). Structures shown as cartoons, colored as in (4c). Intermolecular 
distances to indicate proximity for SARS-CoV PLpro-Ub are shown as black dashes with the 
ones for SARS-CoV PLpro-ChISG15 in purple, (c) Activity of SARS-CoV PLpro mutants 
towards Ub-, and ISG15-AMC. Corresponding data found in Table 1. 

Figure 6 (a) Cartoon representation of the SARS-CoV PLpro complexed with CmISG15 (green). 
Secondary structure of the PLpro is designated, with helices and loops rendered in gold and P- 
sheets rendered in silver. The structural domains of the PLpro are identified, consisting of the 
fingers (white), palm (purple), thumb (green), and Ubl (red) labeled with colored boxes, (b) 
Cartoon representation in wall-eye stereo view, with transparent surfaces, of the interactions 
between the thumb domain (gold) and Ubl-domain (gray) for the SARS-CoV PLpro-CmISG15 
bound structure. Amino acids involved in interactions are represented as sticks and hydrogen 

o 

bonds are indicated by black dashes and distances (A) labeled in red. (c) Ribbon overlay of 
SARS-CoV PLpro complexed with CmISG15 (gold) and SARS-CoV PLpro complexed with 
ChISG15 (reddish orange) demonstrating the movement of the Ubl-domain. (d) Overlay of 
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SARS-CoV PLpro (gold) complexed with CmISG15 (green) and the ChISG15 (blue) from the 
SARS-CoV PLpro complex based on a secondary structure alignment of the respective PLpros. 
The approximate degree of shift in orientation of analogous a-helices is indicated, (e/f) 
Comparison between the SARS-CoV PLpro (reddish orange) complexed with ChISG15 (blue) 

o 

and CmISG15 (green). Electrostatic interactions indicated with black dashes and distances (A) 
labeled in red or yellow, respectively, (g/h) Interactions between SARS-CoV PLpro (gold) and 

o 

CmISG15 (green) with electrostatic interactions indicated with black dashes and distances (A) 
labeled in red and hydrophobic interactions indicated by the yellow transparent surface. 

Figure 7 (a) Secondary structure of monomer A of full-length mISG15 (purple) with P-sheets 
colored in silver, (b) Monomers of mISG15 were overlaid, aligning the C-terminal domains, to 
observe the flexibility of the monomeric units, (c) Monomers of hISG15 (green; PDB 1Z2M, 
gray; PDB 3R66, blue; PDB 3PSE) were overlaid, aligning the C-terminal domains, to observe 
the flexibility of the monomeric units. 

Figure 8 (a) Overlay of mISG15 (purple) and hISG15 (green; PDB 1Z2M) shown in ribbons, (b) 
Comparison of interactions responsible for the structural conformations of hISG15 (green) and 

o 

subunit A mISG15 (purple) with dashed lines representing bond lengths in A. (c) Interactions 
observed when mISG15 is forced into the hISG15 conformation (left panel) and vice versa (right 
panel), (d) Electrostatic surface of mISG15 and hISG15. 

Figure 9 (a) SARS-CoV PLpro (reddish orange) from the ChISG15 bound structure with full- 
length human (green; PDB 1Z2M) overlaid with full-length human based off of a secondary 
alignment using Ub (light gray; PDB 4MM3). The distance between the two-overlaid full-length 

o 

human structures is measured in A and is indicated by a blue dashed line. For all figures the 
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highlighted yellow portion demonstrates the secondary binding site ridgeline, (b) SARS-CoV 
PLpro (gold) from the CmISG15 bound structure with full-length human (green; PDB 1Z2M) 
based off of a secondary structure alignment using ChISG15. (c) Close-up overlay of the 
secondary binding site helix from SARS-CoV PLpro (gold) from the CmISG15 bound structure, 
SARS-CoV PLpro (reddish orange) from the ChISG15 bound structure, and SARS-CoV PLpro 
(dark gray; PDB 5E6J) from the K48 di-Ub bound structure. The distance between the helix 
from SARS-CoV PLpro (gold) from the CmISG15 bound structure and the helix from SARS- 

o 

CoV PLpro (reddish orange) from the ChISG15 bound structure is measured in A and is 
indicated by a blue dashed line, (d) SARS-CoV PLpro (reddish orange) from the ChISG15 
bound structure with full-length human (green; PDB 1Z2M) overlayed with K48 di-Ub (dark 
brown; PDB 5E6J) both based off of a secondary alignment using ChISG15. (e) SARS-CoV 
PLpro (gold) from the CmISG15 bound structure with full-length mouse overlayed with full- 
length human (green; PDB 1Z2M) both based off a secondary structure alignment using 
CmISG15. 
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Table 1 . Kinetic Parameters of SARS PLpro WT and mutants with different ubiquitin-based fluorescent 


substrates. 


Substrate 

SARS PLpro Enzymes 

Kinetic Parameter 

WT 

R167E 

M209A 

Q233E 

RLRGG-AMC 3 





kcat/Km (pM' 1 min' 1 ) 

0.141+0.002 

0.062 + 0.001 

0.727 + 0.007 

0.211+0.002 

Ub-AMC a 





WKmiliM' 1 min' 1 ) 

7.22 + 0.56 

58.0 + 2.5 

2.90 + 0.32 

3.14 + 0.32 

ISG15-AMC b 





kca/KmiliM" 1 min' 1 ) 

50.7 + 9.0 

2.76 + 0.10 a 

54.3 + 12.8 

132 + 23 


J Best-fit slope values derived k app for nonsaturating substrates and approximates kc at /K m 
b Steady-state values following Michaelis-Menten kinetics are reported from duplicate measurements 
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Table 2. Isothermal Titration Calorimetry. 


Protein 

N 

k d 

AH 

AG 

-TAS 


(sites) 

(pM) 

(kJ/mol) 

(kJ/mol) 

(kJ/mol) 

SARS-CoV PLpro 






Human 

0.932 + 0.032 

20.50 + 4.48 

-27.20+ 1.90 

-26.80 

0.38 

C-Human 

1.31+0.015 

57.6 + 3.21 

-4.03 + .019 

-24.18 

-20.17 

Sheep 

1 

200 + 41 

-106 + 10 

-21.10 

84.8 

Mouse 

1 

196+68 

-92 + 15 

-21.20 

70.4 

C-Mouse 

1 

1870+ 173 

-20.71 +0.486 

-15.56 

5.10 

MERS-CoV PLpro 






Human 

0.785 + 0.027 

59.3 + 12.7 

-11.34 + 0.268 

-24.15 

-12.80 

Sheep 

1 

172 + 38.5 

-19.92 + 0.420 

-21.51 

-1.59 

Mouse 

1 

376 + 52.8 

-43.51 +0.535 

-19.54 

24.10 
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Table 3. Data collection and refinement statistics. 



SARS-CoV 

PLpro-ChISG15 (PDB 
entry 5TL6)* 

SARS-CoV 

PLpro-CmISG 15 (PDB 
entry 5TL7)* 

mISG15 (PDB entry 
5TLA)* 

Data collection 




Space group 

P2 1 2 1 2j 

P2 1 2 1 2 1 

P2 

Cell dimensions 




a, b , c (A) 

46.9, 87.0, 221.5 

76.0, 98.2, 106.6 

85.8 60.2 172.7 

a, P,y (°) 

90, 90, 90 

90, 90, 90 

90, 93.1,90 

Resolution (A) 

50.00-2.62 (2.67-2.62) t 

50.00-2.44 (2.53-2.44) t 

50.00-3.25 (3.31-3.25)' f 

^sym Emerge 

11.1 (58.3) 

9.9 (74.2) 

5.2 (37.5) 

//o/ 

12.7 (2.07) 

18.0(2.71) 

17.4 (2.27) 

Completeness (%) 

98.4 (96.0) 

99.1 (99.8) 

98.4 (88.0) 

Redundancy 

5.0 (3.9) 

2.4 (2.1) 

2.4 (2.1) 

Refinement 




Resolution (A) 

38.67-2.62 (2.71-2.62) 

32.20-2.44 (2.53-2.44) 

42.18-3.25 (3.36-3.25) 

No. reflections 

27,644 

29,886 

27,929 

T^work (%)/ /'’free (%) 

19.1/25.0 

18.7/25.3 

25.6/30.9 

No. atoms 




Protein 

6177 

6218 

11,491 

Ligand/ion 

15 

10 

0 

Water 

133 

231 

0 

B-factors 




Protein 

49.27 

61.13 

86.7 

Ligand/ion 

59.94 

57.02 

0 

Water 

41.95 

58.23 

0 

R.m.s. deviations 
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Bond lengths (A) 0.004 0.008 0.002 

Bond angles (°) 0.54 1.1 0.57 

*Datasets collected from a single crystal for each structure. ^Values in parentheses are for highest-resolution shell. 
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Highlights 


• Coronavirus PLPs possess immune modulating deubiquitinating/delSGylating activity 

• Different species’ ISG15s have different binding characteristics with PLPs 

• New structures of PLPs with an ISG15 domain expose different binding orientations 

• A structure of full length mouse ISG15 reveals overall differences from human ISG15 

• ISG15 variation affects PLP function and potentially host immune responses 
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